Magnetic Frustration and Iron- Vacancy Ordering in Iron-Chalcogenide 
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We show that the magnetic and vacancy orders in the 122 (Ai- y Fe2- x Se2) iron-chalcogenides 
can be naturally derived from the J\ — Ji — J3 model with J± being the ferromagnetic (FM) nearest 
neighbor exchange coupling and J2, J3 being the antiferromagnetic (AFM) next and third nearest 
neighbor ones respectively, previously proposed to describe the magnetism in the ll(FeTe/Se) sys- 
tems. In the 11 systems, the magnetic exchange couplings are extremely frustrated in the ordered 
bi-collinear antiferromagnetic state so that the magnetic transition temperature is low. In the 122 
systems, the formation of iron vacancy order reduces the magnetic frustration and significantly in- 
creases the magnetic transition temperature and the ordered magnetic moment. The pattern of the 
245 iron- vacancy order (y/5 x y/h) observed in experiments is correlated to the maximum reduc- 
tion of magnetic frustration. The nature of the iron- vacancy ordering may hence be electronically 
driven. We explore other possible vacancy patterns and magnetic orders associated with them. We 
also calculate the spin wave excitations and their novel features to test our model. 
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The recent discovery of a new family of iron- 
based superconductors, the 122 iron-chalcogenides 
A(K, Cs, Rb)yFe2-xSe2\ui<M, with a superconducting 
transition temperature even higher than 40 K has at- 
tracted many research attentions. The compounds 
are heavily electron doped with only electron Fermi 
pockets which are mainly located at the M point of 
the unfolded Brillouin zone as shown by both angle- 
resolved photoemission spectroscopy (ARPES)[4|-[6j and 
LDA calculations Among all of iron-based su- 

perconductors, the new materials also have the high- 
est magnetic transition temperature around 500K and 
the largest ordered magnetic moment around 3/i^[lQl.[ll|. 
The materials are also featured with intrinsic iron vacan- 
cies which can order themselves below some transition 
temper at uresjl Oi l 12l - ll6| . Currently, it is unclear how the 
magnetic, vacancy and superconducting orders are cor- 
related. 

For iron-pnicitides, their parental compounds dis- 
play an universal collinear antiferromagnetic (CAFM) 
phase [17], which can be described by a simple magnetic 
exchange model including the nearest neighbor (NN) 
J\ (note, in magnetically ordered state, the effective 
exchange coupling becomes J\ a and Jib due to rota- 
tion symmetry breaking) and the next nearest neigh- 
bor (NNN) J2IH Gjjj. However, for iron-chalcogenides 
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different magnetic orders have been observed jiol 
For example, the 11 iron-chalcogenides, FeTei- x Se Q 
which can achieve the superconducting transition tem- 
perature around 40K under pressure [23|, can display 
both commensurate and incommensurate magnetic state. 
The ordered magnetic moment of the parent compound 
FeTe is about 2.0/ib, significantly larger than iron- 
pnictides as well. Recently, it has been shown that the 
magnetic state of FeTe, which is a bi-collinear antiferro- 



magnetic (BAFM) can be described by a strongly frus- 
trated magnetic model including the nearest neighbor 
(NN) Ji, the next nearest neighbor (NNN) J2 , and the 
third nearest neighbor (TNN) J3, i.e. the J\ — J2 — J3 
mo del 0, [Hj. It has been determined that J\ is ferro- 
magnetic (FM) while J2 and J3 are antiferromagnetic 
(AFM) [26]. It is also interesting to note that in the 
BAFM state, the J\ effectively becomes J\ a and Jib due 
to symmetry breaking, which also help to stabilize the 
BAFM phase. The values of the magnetic exchange cou- 
pling parameters suggest that FeTe is close to a bound- 
ary between an incommensurate magnetic phase and the 
BAFM phase so that the later can lose its stability by 
an introduction of a few percentages of additional Fe 
atoms (20). 

However, it is still hotly debated about whether the 
effective models with local magnetic exchange couplings 
are the right models to describe the magnetism in iron- 
based superconductors. For iron-pnictides, there is less 
debate since the model is simple and has gained support 
from different families of iron-pnictides. However, for 
FeTe, the model lacks of independent verification. Since 
the local physics of the 122 iron-chalogenides should be 
similar to the one of FeTe, it is naturally expected that 
both families of iron-chalogenides should be described by 
similar models. 

In this paper, we show that the magnetism in both 
families of iron-chacogenides can be unifiedly understood 
within the previous model obtained for FeTe. The 245 
(i^Fe^es) pattern (y/h x y/h) of the vacancy orders in 
the 122 system are naturally obtained by maximizing the 
magnetic energy saving and reducing the magnetic frus- 
tration in the model. The magnetic order in the 122 sys- 
tems becomes robust and the ordered moment is signif- 
icantly enhanced. The magnetic transition temperature 
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around 500k of K 2 Fe4 : Se 5 quantitatively agrees with the 
pure magnetic energy gain, comparing to the magnetic 
transition temperature around 70K in FeTe where the 
magnetic frustration is strong. Therefore, the vacancy 
ordering may be originated electronically. We explore 
other possible vacancy patterns and magnetic orders as- 
sociated with them. Moreover, we calculate spin wave 
excitations and identify their novel features to quantita- 
tively test the proposed theory. 

Before we start to consider the effect of vacancy, let's 
review the J\ — J 2 — J 3 model, which was used to describe 
the magnetism in FeTe 2J, [25j . The Hamiltonian is given 

by 



Oi <ij>a 



(1) 



where a = 1,2,3. The BAFM state in FeTe, which 
breaks the rotation symmetry, is stabilized with a cou- 
pling to the small lattice distortion. In the BAFM state, 
this broken symmetry can generate the anisotropy among 
J\ and J 2 couplings. For FeTe, by fitting the spin wave 
spectrum, we have obtained the values of the magnetic 
exchange couplings with J\S ~ — 34mev, J2S ~ 22mev 
and J3S ~ 7mev[26], where S is the spin of each site. 
The anisotropy is shown mainly in Ji, with Jla ~ Jlb ~ 
— 16mev^2(^. The most important feature revealed from 
the fitting is that the NN coupling J\ is FM while the 
other two, J 2 and J3, are AFM. The sum of these three 
coupling parameters is close to zero, which indicates the 
model is extremely magnetically frustrated. 

Let's first ignore the anisotropy of the magnetic ex- 
change coupling Ji caused by symmetry breaking in the 
BAFM state. We focus on the pure magnetic model 
in the tetragonal square lattice. In the BAFM state as 
shown in figHJ both J\ and J2 are frustrated couplings 
and the magnetic energy is only saved by J3. The saved 
magnetic energy per site is given by Ebafm = —2J^S 2 . 
Now, we ask a simple question, what are the possible 
stable magnetic orders if we can remove spins and form 
vacancy patterns in the above model within the similar 
parameter regions? (note, in the following we will take 
S = 1 for convenience.) 

In order to create a vacancy pattern that minimizes 
the magnetic exchange energy, the pattern must suffi- 
ciently suppress magnetic frustration and maximize en- 
ergy saving from each magnetic exchange coupling si- 
multaneously. It is straightforward to show that such 
a pattern is exactly the 245 pattern ( named accord- 
ing to K2Fe^As^) as shown in figEJ which has been 
identified by neutron scattering experiments [To| . In the 
245 vacancy orders, each site has three NN and NNN 
couplings to its neighbors. Moreover one spin in each 
connected five spins (5-1 pattern) along Fe-Fe direction 
is removed. The odd number of links for both NN 
and NN couplings reduce the magnetic frustration com- 
pletely. The 5-1 pattern saves energy from J 3 as well. 
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FIG. 1: The bi-collinear antiferromagnetic phase observed in 
the 11 (FeTe/Se) systems. 



Therefore, the 245 state maximizes the energy saving 
from all magnetic exchange couplings. The magnetic en- 
ergy saving per spin F245 = (— J\ + J 2 + 2Js)/2(we set 
J3 = between two NNN sites if there is vacancy between 
them) . If we ignore the lattice distortion and use the ex- 
change parameters derived from FeTe, F245 ~ 70mev, 
comparing to Ebafm ~ 14mev. Considering the fact 
that the ordered magnetic moment is also around 1.5 
times larger in K2F e^As^ than in FeTe, the ratio of the 
magnetic energy savings between two states is given by 
1.5 x F245/ Ebafm ~ 7.5. This ratio is in a quantita- 
tive agreement with the ratio of their magnetic transi- 
tion temperatures 500k /70k ~ 7.1. (Note: in general, 
the magnetic transition temperature in a quasi two di- 
mensional material is given by T/v oc t^tt, where J c is 

J c 

the coupling strength between layers and E is the in- 
plane magnetic energy saving.) The vacancy pattern is 
hence a result of the combination of energy saving from 
three magnetic couplings including the large contribu- 
tion from the NN ferromagnetic exchange coupling. If 
we consider the symmetry breaking of the state which 
can modify the magnetic coupling strength as shown in 
figHJ the total energy saving in general can be written as 



^245 = (-2 Ji -J[ + 2^2 - J 2 + 2 J 3 - J3V2 where the 



parameters are defined in figEJ 

It is also interesting to explore other possible vacancy 
patterns and their associated magnetic orders. For such 
an frustrated J\ — J2 — J3 magnetic model, there is no 
other pattern, except the 245 ordering, which can get rid 
of the magnetic frustration completely. With the values 
of the exchange couplings given above, we list some pos- 
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FIG. 3: (a) The armchair dimer crystal pattern. (b)The 
square dimer crystal pattern. 

sible patterns (we again set J3 = between two NNN 
sites if there is vacancy between them): (1) Armchair 
dimer crystal patterns (the 212 pattern if we consider 
K 2 FeSe2) shown in ngj3ja). In this pattern, half of spins 
are replaced by vacancies. The magnetic energy saving 
per site is given by Eadc = J 2 — Ji/2. The pattern saves 
both energy from the NN FM and NNN AFM couplings; 
(2) Square dimer crystal patterns shown in ngU^b). The 
energy saving from this structure stems from the NN FM 
Ji and the NNN AFM J 2 , E S dc = ~ 4Jl 5 +2J2 ; (3) The 
234 vacancy patterns (named as ^Fe^Se^) as shown 
in figlU There are two different patterns. The first 
pattern saves energy from J\ and J2 but pays energy 
cost for J3. The energy saving per site is given by 
E234a = — Ji/3 + 2J2 — 2J3. This magnetic structure is 




FIG. 4: The two different possible patterns (234) for 
Fei. 5 Se2. 

similar to the (0, 7r) collinear AFM (CAFM) observed in 
the parental compounds of iron-pnictides. In the second 
pattern, the magnetic state is also similar to CAFM. The 
energy saving is given by E 2 34b = §(^2 — J3/2). In gen- 
eral, the 234 phase where the iron concentration is close 
to 1.5 is most likely supporting a CAFM magnetic phase, 
similar to iron pnictides. In both patterns, there are fer- 
romagnetic moments for each layer. Between layers, they 
are antiferromagnetically coupled. The 234 vacancy pat- 
terns have been observed experimentally [28]. However, 
the magnetic order has not been identified. (4) The 446 
pattern (K^Fe^Seo) as shown in fig|5l There are also two 
possible patterns. The first pattern is a dimer vacancy 
ordering. The magnetic pattern is also close to CAFM 
phase with an energy saving given by £?234a = j2 ~ J |~ 3Jl . 
The second pattern can be viewed as a diagonal stripe 
similar to undoped cuprates[29]. For such a stripe pat- 
tern, there is no J3 coupling. The magnetic structure 
can be predicted to be incommensurate along the stripe 
direction and AFM along the direction perpendicular to 
the stripe. The total energy saving per site is given by 
£?446 = I J 2 + 3j^- The spin angle between two NN sites 
along the stripe direction is given by cosO = — jj- . 

Properties of the 245 state in the large S limit: The 
245 state is clearly the most stable magnetic state in our 
model. In the following, we focus on this state. First, 
we analyze the spin excitation in the large S limit. We 
consider a general model as shown in Fig|2] To start a 
spin wave calculation, we denote the spin sites as follows. 
A generic position of the spin is given by r = mli + 71I2 + 
di, where m, n are integers and li = (2x-y)/V / 5,l 2 = 
(x+2y)/v / 5, di = 0, d 2 = x, d 3 x+y, d 4 y, where 
x, y are unit vectors in the original tetragonal lattice. 
We take the Holstein-Primakoff transform for the given 
ground state: 
For rn + n =even: 

5+(r) = V2S ai (R), (2) 
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FIG. 5: The two different possible patterns for Fex.^Se^'. 
(a) The dimer vacancy ordering similar to CAFM magnetic 
phase, (b) The diagonal stripe pattern where the magnetic 
ordering is incommensurate. 



S-(r) = V2Saf(R), 
S z (r) = S-at(R)oi(R); 



For rn + n =odd: 



5+(r) = 
S-(r) = 
S z (r) = 



25at(R), 
2Sa;(R), 
5 + ot(R)oi(R). 



(3) 



The full Hamiltonian of the model can be 
put in a matrix form. Define ip^(k) = 

(a\(k),al(k),al(k),al(k),ai(-k),a 2 (-k),as(-k),a4 : (-k)), 
and we have 



H 



A(k) B(k) 
B(k) A{k) 



(4) 



J 



A(k) and B(k) are four-by-four matrices, defined by: 



/ Eq J\e lkx J 2 e lkx+lky + J' 3 e l2kx 



A(k) = S 



J\e lky 



E 



Jie iky 

Eq 



f £2 c 2i%\\j / u 



Jke 



\ ■ 



Jie lkx 
Eo 



(5) 



/ J' e - ikx+ik y 



B{k) = S 



J e 2iky 



j/ e ik x ik. 



J' x e~ lk y j^ e -tk x -tk y j se 2ik x \ 

Zy + Jse 2ikx J[e ikx 

J^ e ik x -ik y J 3e 2iky 

/ 



(6) 



where E = -2J X - J 2 + J[ + 2J' 2 + 2J 3 - J3. The lower 
triangle elements are suppressed because both matrices 
are hermit ian. 

By diagonalizing this Hamiltonian, we have 

H= ]T hl(kH(k) + l/2)cui(k), (7) 

i=l,2,3,4;/c 



and 



Oi(k) = ^2u ij (k) lj (k) + Vij (k)^(-k). 



(8) 



where U and V are 4x4 matrices to diagonalize the 
Hamiltonian. The spin reduction due to quantum fluc- 
tuations on each site is given by: 

fa= ^2E/ o dk *J o dk y\Vij(k x ,k y )\ 2 . (9) 

It is also interesting to calculate the dynamic factor 
S(u) = — 9[x(cj)], wnere the local susceptibility x(cj) 



is defined as 



E 

q 



J c/re-"-[(5 + ( 9 ,r)5-(- 9 ,0)) 



+<£-(«, t)S + (- 9) 0))]/JV. 
After taking Fourier transform, we obtain 



(10) 



k i,j 

- Ej(k) + i0 + Lj + Ej(k)+iO^ N ' ^ 

Due to the enlarged unit cell, the spinwave BZ (SBZ) 
is smaller than the original unfolded BZ for one Fe 
per unit cell. The SBZ is cornered by the following 
four (equivalent) points (— 37r/5, — 7r/5), (7r/5, — 37r/5), 
(37r/5, 7r/5), and (— 7r/5, 37r/5), and these points are 
called M s points, different from the real M point in 
lattice BZ. It is convenient to define X s = (27r/5, — 7r/5) 
and Y s = (— 7r/5, — 2tt/5) for future discussion. The 
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FIG. 6: The spinwave spectra of the lowest three 
branches (the other one is too high to be drawn in the 
same plot) using the parameters ( Ji, J{, J2, J2, J3, J3) = 
(-30, -10,20,20,9,0)meV in (a) and ( J u J[ , J 2 , J2, J 3 , J3) = 
(-30, -30, 20, 20, 9, 0)meV in (b). The corresponding band 
structures are given in (b) and (d), in which the read dashed 
lines represent the highest branch and black solid lines the 
lowest three branches. M s , X s and Y s are high symmetry 
points in the spinwave BZ that are defined in text. 

spinwave has four branches. Depending on the param- 
eters, there can be finite gaps between the first and 
the second branches and between the third and the 
fourth ones . The middle two branches are degenerate 
at T and M s points and are in general close to each 
other. At high symmetry points, the spinwave energy 
has analytic expressions. We take S = 1 for conve- 
nie nce. At M s and T: gi(T) = 0,E 2 (T) = E 3 (T) = 
2V(J e -J1- J2- J'M ~ Ji ~ J2 + Je ~ Jj, and 
E 4 (T) = 2V / 2 v / ( J i - 2 Ji)(Je ~ Ji) where J e = J' 2 + J 3 
are the effective antiferromagnetic coupling strength 
between two 4-spin plaquettes. At X 3 , E\{X S ) = 

2(Jl + V( J e -J3- Wl +Je~J'z- Jl)), E 2 (X S ) = 

2yfjf + ( J 2 - 2j\){ J 2 - 2J 3 ) + uj, where u 
J[(J e ~ J 2 ) - + 2J e - 2J 2 ), and £ 4 (X 5 ) = 

2(-Ji + y/(J e - J3 - Ji)(J{ + J e ~ J3 ~ Ji)), which 
also defines the spinwave band width. The first spinwave 
branch effectively describes the spin fluctuations of 
4-spin plaquettes. 

We calculate spinwave spectra and related features for 
three different parameter sets. We first adopt the pa- 
rameters: Ji = (Ji a + Jib)/ 2 = — 30meV, J 2 = 20meV, 
J3 = 9meV. Due to lattice distortion that draws the four 
spins in each new unit cell closer, the primed param- 
eters can be reduced from unprimed ones by the pres- 
ence of iron vacancies: J[ = — lOmeV, J 2 = 20meV and 
J3 = OmeV. For these parameters the spinwave spectra 
have a finite gap Ai 2 = 10.2meV between the first and 
the second branches, and a larger gap A34 = 86.2meV 
between the third and fourth branches. The full band- 



width is W = 167.5meV. The first branch starts from 
zero energy at T and M s and reaches maximum at X s ; 
the second and third branches start from minimum at 
(37r/10, 7r/10) and its three C4 equivalents, and reach 
maximum at (0.457T, 0.127r). The fourth branch starts 
from Ki(r) at F and M s to E^Xg) at X s . The ordered 
moment correction from spinwave is 5m = 0.095. See 
FigE^a,b) for details. 

In choosing the second parameter set, we consider a 
purely magnetic model without lattice distortion. There- 
fore we take J[ = Ji, J 2 = J 2 . However one still has 
Js ^> J3 ~ because of the vacancy along the exchange 
path for J3. For this parameter set, there is no gap be- 
tween the first and second branches, but the gap between 
third and fourth branch remains finite. The full band- 
width in this case is W — 142.7meV. The first three 
branches start from zero energy at T and M s , and reach 
maximum at (37r/10, n/lO) and its three C4 equivalents. 
The fourth branch's minimum and maximum appear at 
the same points as in the first parameter set. The or- 
dered moment correction from spinwave is 5m — 0.197. 
See Fig|6fc,d) for details. 

In FeTe, we found that the third neighbor exchange 
Js is needed to explain both low and high energy re- 
gions of dynamic spin susceptibility observed in neu- 
tron experiment, while a parameter fit without J3 sig- 
nificantly underestimates the spin wave bandwidth. J3 
is also important in the current system to stabilize the 
245 state. For comparison, we adopt the parameters 
we used to fit FeTe system without J3 (see supple- 
mentary materials in Ref . [26] ) : ( J\ , J[ , J 2 , J 2 , J3 , J3 ) = 
(—12, —4, 9, 9, 0, 0)meV. The spinwave correction to mo- 
ment in this case is 5m = 0.43, indicating much stronger 
spin fluctuation. 

The dynamic factor is shown in FigJTl Note that 
around zero energy the dynamic factor should drop to 
zero in a 3D system, but converge to a finite value in our 
2D system. If we include a small J z in the beginning, we 
can make the curve drop to zero at zero energy, however 
the dynamic factor at uj ^> J z is unaffected. See FigJTl 
for the dynamic factor with all three parameter sets given 
above. 

It is also interesting to check what is the closest mag- 
netic state if the exchange coupling parameters can be 
varied in the 245 state. In the reasonable parameter 
space, the following states can be considered: 1. The 
ferromagnetic state if J\ and J[ are large. 2. the CAFM 
state when J 2 is large. 3. The incommensurate phase. 
The ground state energy for the first two are easy to cal- 
culate: E FM = 4 Ji + 2 J[ + 2 J 2 + 4 J' 2 + 4 J 3 + 24 and 
Eafm2 — — 2 J 2 — 4J 2 + 4J3 + 2 J3. Comparing the ener- 
gies, we have if | J[\ > 2(J 2 + J3), the FM phase wins over 
the 245 magnetic phase; and if 2| Ji| — | J[\ < 2J 2 — 4J3, 
CAFM state has lower energy than the 245 magnetic 
state and thus becomes preferable. Most importantly, an 
ICM state is the closest competitor to the 245 magnetic 
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U'l|/J2 



FIG. 7: The dynamic factor S(cj) contributed by spin- 
wave. Black solid line corresponds to (Ji, J{, J 2 , J 2 , J3, J3) = 
(— 30, — 10, 20, 20, 9, 0)meV; blue dashed line corresponds to 
(Ji, Ji, J 2 , J 2 ,Js, J3) = (-30, -30,20,20,9,0)meV; and red 
dotted line (Ji, J{ , J 2 , J 2 , J 3 , J3) = (-12, -4, 9, 9, 0, 0)meV. 
The main panel only includes contribution from lowest three 
branches as the highest branch is too high to be plotted in a 
same energy range. The inset shows the contribution by the 
fourth and highest branch. 



phase. An analysis of the spinwave Hamiltonian around 
T point gives us the phase boundary between the 245 
magnetic state and ICM as drawn in fig|HJ The bound- 
ary is set by the equation: J\{J[ + 2(J e — Jg)) + J[(J2 — 



J e + J^)+2(J 2 (J e 



2J 3 (- + J£) = 0. We draw 



fig |8] (left) for a phase diagram against J\ and J{, fixing 
J 2 = J 2 and J3 = 9J2/2O, and fig|8] (right) for a phase 
diagram against J\ and J2, fixing J[ = Ji, J 2 = J2 and 

4 = 0. 

Discussion: We have shown that the 245 vacancy 
pattern can be naturally obtained through saving mag- 
netic energy in a strongly frustrated magnetic model 
J\ — J2 — J3 model. A few remarks regarding our the- 
ory need to be addressed. First, in our above analysis, 
we only consider the magnetic energy and ignore other 
possible energy sources. In real materials, it is possible 
that other energy sources may make a significant con- 
tribution as well. For example, creating a dimer of two 
vacancies may cost higher or lower energy than being 
two individual vacancies. Therefore, comparing a dimer 
vacancy pattern to a non-dimer pattern requires further 
careful consideration. The LDA calculation may help 
to address it [30]. Second, we expect a lattice distor- 
tion to naturally take place above or equal to the mag- 
netic transition temperature, similar to what has been 
observed in iron-pnictides[17]. Third, some LDA calcu- 
lations have predicted the values of magnetic exchange 
couplings for K 2 F e^As^]^. However, their values are 
not consistent with our picture. In their results, the NN 
exchange coupling J[ is antiferromagnetic and there is 
no J3. We expect that both J\ and J[ should be FM 
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FIG. 8: The phase diagrams against J\ and J[ (left) and 
Ji and J 2 (right). Solid lines mark second order transitions 
calculated from spinwave and dotted lines mark possible first 
order transitions. 



and J3 should be significant. Third, it is very important 
to notice that the NN exchange coupling is strongly fer- 
romagnetic, which contributes the largest energy in the 
245 phase. It leads to a high magnetic transition tem- 
perature observed in the 245 phase. It is interesting to 
note that the 245 pattern has also been obtained in one 
band t — J — V model [3l| which is driven by magnetic 
energy saving in proper parameter regions. Fourth, the 
ferromagnetic coupling does not contribute to any spin 
singlet pairing in superconducting (SC) state. Therefore, 
although the magnetic transition temperature is high in 
these new materials, we do not expect the SC transi- 
tion temperature can be scaled as the magnetic transi- 
tion temperature here because the SC transition temper- 
ature is determined by the strength of the AFM cou- 
plings. Fifth, the SC pairing is mainly determined by 
Jo [32] . We expect the J 2 value is similar to the one in 
iron-pnictieds[33|, |34j since the SC transition tempera- 
tures for both materials are similar. Sixth, the existence 
of J3 distinguishes iron-chacogenides from iron-pnictides 
and generate many different magnetic states. The effect 
of J3 in superconducting state will be addressed in a dif- 
ferent paper. In fact, J3 can enhance the s-wave pairing 
significantly when electron pockets dominate over hole 
pockets. Finally, the 245 pattern is so strong in mag- 
netism. We do not think this state can be coexisted with 
superconducting phase. 

In summary, the J\ — J 2 — J3 model provides a natural 
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and unified understanding for magnetism and vacancy 
ordering in iron chacogenides. The 245 vacancy pattern 
reduces the magnetic frustration and achieves the maxi- 
mum magnetic energy saving. The reduction of frustra- 
tion increases the ordered magnetic moment as well as 
strongly enhancing the magnetic transition temperature. 
(Note, after finishing our paper, we find similar spin- 
wave excitations for the 245 pattern are calculated and 
discussed in[35[ although the parameters in their model 
and the physical motivation are very different from us.) 

Acknowledgement: JP thanks S. A. Kivelson, Hong 
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